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. Weakly coupled superconducting layers are described by the three- 

o , 

dimensional XY model with strong coupling in two directions and weak cou- 

pling in the third direction. For the usual Josephson-type interplane coupling 

the coherence between the layers is lost at the same temperature as that within 

C , the layers. Thus a low-temperature layer decoupling due to a proliferation of 

fluxons between planes, as proposed by Friedel, does not occur in this case. 

^ , However, for a modified interplane coupling there are two phase transitions, 

one of a Kosterlitz-Thouless type from a disordered high-temperature phase 

to an intermediate phase with phase coherence only parallel to the layers, 

the second from this effectively two-dimensional phase to a three-dimensional 

phase with coherence in all directions and a finite "n-state" order parameter. 

Thus we do find a "Friedel transition" for this special class of models. 
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The Berezinskii-Kosterlitz-Thouless (BKT) transition!!!, mediated by the unbinding of 
vortex pairs, has been clearly observed in superfluid films!. More recently, nonlinear trans- 
port experiments in layered high-temperature superconductors^ have also shown typical 
signatures of vortex unbinding slightly below the critical temperature T c . This is surpris- 
ing, since the Josephson coupling between the layers renders the system three-dimensional, 
in particular close to T c . Specific heat experiments for YBCO single crystals indeed give 
evidence for three-dimensional critical behavioiH. 

Some time ago Friedelil has argued that the interlayer coupling could be effectively sup- 
pressed by a proliferation of Josephson vortex loops or ffuxons between the layers. A simple 
estimate for the fluxon energy suggests that the layers are decoupled at a temperature 
T* < T c , so that there would be a two-dimensional regime for T* < T < T c , with a BKT 
transition at T c . Unfortunately, a closer look at the proble ml indicates that this "Friedel 
transition" does not occur below T c i. Therefore the question arises why nevertheless in 
nonlinear transport experiments BKT signatures are clearly observed. A way out of this 
dilemma has been offered by Jensen and Minnhageni, who realized that the Lorentz force 
acting on the vortices can overcome the interlayer confinement of vortex pairs. 

In this Letter we study again the possible existence of a two-dimensional regime below the 
critical temperature, starting from the classical XY model with strong intralayer coupling J\\ 
and weak interlayer coupling J±. A simple criterion for the decoupling transition together 
with Monte Carlo simulations shows that a layer decoupling below T c does not occur, in 
agreement with previous studies!. We attribute this negative result to a strong increase of 
the fluxon energy as a function of temperature. We turn then to the question whether a 
decoupling can be excluded on general grounds, by considering two modified XY models. 
For the first model, representing a superlattice of high and low T c layersS, a decoupling 
seems to occur around the low T c , but a closer look shows that this apparent transition 
is in reality a crossover from strong to weak interplane coherence. In the second model a 
modified interlayer coupling is considered, allowing for n equivalent phase differences for the 
superconducting order parameters of adjacent layers0. In this case a decoupling at T* < T c 



is clearly observed for n > 2. 

Incidentally, the problem is of more general relevance, as the question of long-range co- 
herence arises also in other quantum systems. One can for instance ask whether a disordered 
layered system of electrons, with a large difference between the masses for the motion parallel 
and perpendicular to the layers, can be metallic in one and insulating in another direction. 
Anderson has argued that this can happen, if electron-electron interactions are taken into 
account@. A measure for quantum coherence in the case of electronic transport is the Drude 
weight or charge stiffness, while in the context of superfluidity the relevant quantity is the 
superfluid density p s (or the helicity modulus, in the language of the XY model). 

We consider the classical XY model on a cubic lattice 

H = -J2 J » C0S (Vi-<Pi+»), (!) 

i,fj, 

where \i = x,y,z, J x = J y = J\\, J z = J±, and the phases are restricted to < <fi < 2n. 
This model can be derived from the anisotropic Ginzburg-Landau or the Lawrence- Doniach 
model by neglecting both the fluctuations of the electromagnetic field and the amplitude 
fluctuations of the order parameter. There are good arguments for doing this , although 
these degrees of freedom may become relevant within the critical regionllllii. The helicity 
modulus is defined as the second derivative of the free energy with respect to a constant 
phase gradient in the direction /i, and can be written as 

According to Friedel's original suggestion!, for very weak interlayer coupling Yj_ would vanish 
at a lower temperature than Tii, leaving an intermediate temperature region of essentially 
2d character. 

We present now a simple argument against a layer decoupling below the critical temper- 
ature, by expanding Tj_ in powers of J±. The leading order coefficient (~ J]_) turns out to 
be the difference of two equal contributions, each of them given by S = / d 2 rc 2 (r). The 2d 
correlation function c(r) = (cos(<£>j — y?i+r)) decays exponentially as a function of distance r 
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for T > Tkt- Therefore the leading term in the expansion of vanishes identically above 
the BKT transition (and the same can be shown for all higher order terms). For T < T^t 
c(r) ~ r~ r *( T \ where rj(T) ~ 1/4 00, which implies that the quantity S is infinite below 
the BKT transition. We tentatively associate the temperature where S diverges with the 
transition from an effectively 2d phase of decoupled layers to a 3d phase with finite Tj_. Ac- 
cording to this criterion the decoupling transition temperature T* coincides with the critical 
temperature T c = Tkt in the limit J± — > 0. 

In order to gain a qualitative understanding of the temperature dependent helicity mod- 
ulus, we first use the renormalized harmonic approximation (RHA), where the Hamiltonian 
(H) is replaced by an effective harmonic term 

H= -YsJ^Vi-Vi+vY (3) 

with variational parameters J M . These parameters turn out to be identical to the helic- 
ity moduli, J M = T M . The anisotropy Tj_/T|| diminishes with increasing temperature, but 
remains finite up to the critical temperature, where both Yj_ and Tii drop to zero discon- 
tinuously (Fig.l). We have also performed Monte Carlo simulations using the standard 
Metropolis algorithm. The results presented in Fig.l show that the jumps obtained in the 
RHA are artifacts and that the helicity moduli tend to zero continuously. We note the excel- 
lent agreement between the two methods at low temperatures, as expected. Our numerical 
results for J± = 0.1 Jy are consistent with a simultaneous loss of coherence parallel and 
perpendicular to the layers at T c « 1.33 Jn. Therefore the temperatures T* and T c coincide, 
in agreement with the simple criterion discussed above. 

We focus now our attention on the role of fluxons by considering an approximate version 
of the XY model, where the nonlinearity is retained only for the interlayer coupling, 

H = ~j (fi- Vi+v) 2 - J± cos (^ - Vi+z) ■ (4) 

This is an excellent approximation for the original XY model for T Jn. The helicity 
modulus parallel to the layers is constant and given by Ty = Jn. In order to calculate Tj_ 
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we treat the nonlinear term using the Villain approximation^, which is very accurate both 
at low (T JjJ and at high temperatures (T 3> J^)lll. The partition function is factorized 
into a term representing the in-plane harmonic fluctuations and the "fluxon contribution" 

Zfl = e - 2 * 2f3 ^ miVijmj , (5) 

{mi} 

where the variables mi are integers. The interaction V%j is the Fourier transform of 

v ,s = j* 2 - cos q x - cos q y 

± 2 — cosq x — cosq y + (J±/ J\\)(l — cosq z ) ' 

where the effective interlayer coupling is given by 

Ji= ( 2/3i ° g ^r <7) 

for /3Jj_ <C 1. We notice that Vij is exactly equal to the interaction energy of two elementary 
fluxons calculated in the usual vortex loop representation of the 3d XY modelEl. The 
variables are the quantum numbers of the fluxons, and large loops can be constructed 
by adding elementary fluxons. Since the energy scale J]_ increases roughly linearly with 
temperature the multiplication of fluxons is strongly slowed down. The helicity modulus 
perpendicular to the layers, Tj_, can be expressed in terms of fluxon variables, 

-N 



T ± = Jl<|l-47r 2 /3Jli-(fE^] )}■ (8) 



This confirms that for this approximate model the proliferation of fluxons is directly related 
to the loss of interlayer coherence. Eq. (|j) could in principle be used for calculating Tj_, but 
since the couplings decrease only slowly with distance (like a dipole-dipole interaction) 
we have calculated starting from the original expression ([|), using both the RHA and 
Monte Carlo simulations. The results, shown in Fig. 2, demonstrate that even for small 
couplings Jj_ the helicity modulus vanishes only far above the BKT transition, i.e. in 
a temperature region where the Hamiltonian (|]) is no longer a good approximation for the 
original XY model. Nevertheless, the model defined by Eq. (Q) is interesting, since it does 
show a layer decoupling without simultaneous loss of intraplane coherence. The transition 



temperature T* for infinitesimal interplane coupling Jj_ can be easily calculated from the 
perturbation expansion described above. In the present case the correlation function c(r) 
decays algebraically with an exponent i](T) = T/(2irJn), and T* is simply given by the 
relation rj(T*) = 1. Fig.2 shows that the resulting value T* = 2tcJ\\ is consistent with the 
Monte Carlo data. The vanishing of Tj_ for T > T* implies that also (cos(<^)) vanishes, 
although the susceptibility remains infinite up to 2T*. 

The analysis given above strongly supports the view! that the interlayer coherence is not 
destroyed by thermal fluctuations below the critical temperature in the 3d XY model, even 
for arbitrarily small interlayer coupling. We now address the question whether this conclu- 
sion is generally valid, by studying two modifications of the original XY model. The first 
represents a stack of layers with periodically varying mtralayer couplings. The extensively 
studied superlattices of low and high T c layers are nice realizations of such a mode!0. The 
second modification concerns the interlayer coupling, which is replaced by a different form, 
namely J± cos(y?j — (pi +z ) is substituted by (J±/n 2 ) cos[n(^ — (fi+ z )), where n is an arbitrary 
positive integer (> 2). Although such a term cannot be excluded within Ginzburg-Landau 
theory, it would in general be expected to be much smaller than the conventional term with 
n = l. 

We consider first a superlattice consisting of one " strong" layer with intralayer coupling 
Ji| , alternating with n "weak" layers with < Jy , and a constant interlayer coupling 
J_i_. The Monte Carlo results (Fig. 3) show that nothing spectacular happens for n = 1. 
For n = 3, however, the helicity modulus Tii exhibits a kink slightly above the low T c (f» 
J|j 2 ' ) ), while simultaneously drops practically to zero. There is apparently a region 
with vanishing interlayer coherence between this temperature and the critical temperature 
T c w Jy 1 ^. However, a true decoupling is very unlikely. In the extreme situation where 

J \\ = °' 

one can integrate out the variables of the weak layers and deduce a model involving 
only the strong layers with an effective interlayer coupling J±* w (/3J±/2) n J±. We can then 
use our previous results for one type of layers and conclude that the helicity modulus 
remains finite (though very small) up to T c . Thus a Friedel transition does not occur for 
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this type of superlattices. 

We turn now to the second modification where the interlayer interaction is replaced by 
(J±/n 2 ) cos[n (ifi — (pi +x )]. Monte Carlo results for n = 2, 3, 4 are shown in Fig. 4. For n > 2 
the helicity modulus Yj_ drops to zero far below the critical temperature where Ty vanishes, 
while for n = 2 both transitions seem to occur at the same temperature. This observation 
is in agreement with the analysis of perturbation theory. In the present situation we have 
to consider the correlation function 

c n (r) = (cos[n - Vi+r )]) ~ r"^ (9) 

For the exponent r] n we use the relational rj n = n 2 i] together with the numerical values for 
?7(T)0. The criterion t] n (T*) = 1 yields decoupling temperatures T* in good agreement with 
those found in the Monte Carlo simulations for n = 3,4. For n = 2 the reported value 
t](Tkt) = 1/4 indicates that both transitions - the Friedel and the BKT transition - occur 
simultaneously. 

In summary, our numerical simulations and perturbative arguments confirm that for the 
layered XY model an intermediate effectively 2d phase does not exist, even for arbitrarily 
small interlayer couplings, J± <C J\\. This absence of a low-temperature decoupling transition 
is nicely illustrated in a simplified version of the model where the role of fluxons is particularly 
transparent. The interplane coherence also persists for superlattices of high- and low-T c 
layers, although an apparent decoupling is observed for thick enough low-T c layers. In 
contrast, a low-temperature decoupling transition is found for an interlayer coupling with 
n-fold symmetry, at least for n > 2. 

We gratefully acknowledge useful discussions with P. W. Anderson, L. N. Bulaevski, D. 
Feinberg, B. Horovitz, S. E. Korshunov and S. Shenoy. This work has been supported by 
the Swiss National Foundation through grant Nos. 4030-32799 and 20-40672.94. 
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FIGURES 

FIG. 1. Helicity moduli Ty (full circles) and (open circles) of the anisotropic XY model 
with J±/J\\ = 0.1 as functions of temperature as obtained from Monte Carlo simulations of a 
36 x 36 x 36 lattice. The solid curves show the RHA results. 

FIG. 2. Helicity modulus of the approximate version (Eq. (4)) of the XY model with 
J±/J\\ = 0-1 (open circles) and J±/J\\ = 0.5 (full circles). The symbols are Monte Carlo data and 
the full lines are RHA results. Arrows indicate the Kosterlitz-Thouless temperature Tkt and the 
layer decoupling temperature T* = 2nJ\\, respectively. 

FIG. 3. Helicity moduli Ty (full symbols) and (open symbols) for the superlattice model 
(as defined in the text) with / J\\ ] = 0.3 and J±/J^ = 0.1. Monte Carlo data for a 36 x 36 x 36 
lattice with one (circles) and three (triangles) weak layers are presented. The inset shows Tj_ in 
the crossover region. 

FIG. 4. Helicity moduli Ty (full symbols) and (open symbols) of the XY model with 
modified interplane interaction (Jj_/n 2 ) cos[n (ifi — ipi+ z )] with J±/J\\ = 0.1. Monte Carlo data for 
n = 2 (circles), n = 3 (triangles) and n = 4 (diamonds) are shown. The intersection of the dashed 
line with slope 2/tt and the Ty curves locates the BKT transition at Tkt- Arrows indicate the 
calculated layer decoupling temperatures T* for n = 3,4. 
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Fig. 2 
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